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The  p-version  of  the  finite  element  method 
/  for  constraint  boundary  conditions 
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Abstract .  The-  paper  addresses  the  Implementation  of  general 
constraint  boundary  conditions  for  a  system  of  equations  by  the 
p-verslon  of  the  finite  element  method.  By  constraint  boundary 
condi ti onsN*eT mean  conditions  where  some  relation  between  the 
components  Is  prescribed  at  the  boundary.  Optimal  error  bounds 
are  proven.  t  On j  5 x 
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There  Is  a  large  variety  of  boundary  conditions  for  systems 


of  differential  equations  of  elliptic  type.  Some  physically  natu¬ 
ral  conditions  may  be  formulated  by  a  variational  approach  through 
constraint  conditions.  For  example,  the  two  dimensional  elastici¬ 
ty  problem  can  be  formulated  as  the  minimization  of  a  quadratic 
functional  F(u) ,  u  •  (Uj,u2)  over  a  set  H  satisfying 

(Hj(Q))2  c  H  c  (H* (Q ) ) 2 . 

’Selections  of  H  then  characterize  the  boundary  conditions. 

Obviously  the  choice  H  ■  (Hg(Q))2  induces  the  (essential) 

Dirlchlet  conditions,  l.e.,  the  displacement  is  given  on  dQ , 

1  2 

while  H  *  (H  (Q ) )  induces  the  (natural)  Neumann  conditions, 
l.e.,  the  tractions  are  prescribed  on  dQ .  In  addition  to  these 
classical  conditions  other  types  are  important  in  applications. 

One  of  these  conditions  is  characterized  by 

(1.1)  H  -  ( ( Uj , u2 )  c  (H1(Q))2|u1^1(s)  ♦  u2*2(s)  -  0  on  dQ ) 

where  end  *2  0iven  functions  defined  on  dQ .  These  con- 

dltlone  are  in  the  most  simple  case  the  symmetry  conditions  and  in 
general  traction  free  constraints  at  the  boundary. 

So  far  we  have  only  mentioned  homogeneous  boundary  conditions. 
Nonhomogeneous  conditions  are  defined  in  the  usual  way,  when  the 
minimization  of  F  is  over  a  hyperplane  Kvs  (u  +  v|  u  e  h, 
v  c  (H* (Q ) ) 2 ) . 

The  constraint  boundary  condition  we  mentioned  above  is  a  type 
of  essential  condition.  Hence  when  solving  such  problems  by  the 
finite  element  method  in  general  and  by  the  p  or  h-p  versions 


In  particular,  we  face  the  problea  of  implementing  the  nonhomogen- 
neoue  boundary  conditions  (which  are  outside  the  finite  element 
space ) . 

The  p  and  h-p  versions  are  recent  developments,  where  p, 
the  degree  of  the  elements  used  is  not  fixed  but  is  increasing. 

This  Is  In  contrast  to  the  classical  h-verslon,  where  the  degree 
p  is  kept  fixed.  The  first  commercial  programs  available  are 
PROBE  (Noetic  Tech.,  St.  Louis)  and  FIESTA  (ISMES,  Bergamo,  Italy). 

The  implementation  of  Dirichlet  boundary  conditions  for  the 
p-version  of  the  finite  element  method  has  been  addressed  by  us  in 
[2]  and  [4].  A  general  survey  on  the  state  of  the  art  of  the  p 
and  h-p  versions  nay  be  found  in  [1]. 

In  this  paper  we  will  address  the  implementation  of  the  con- 

o 

stralnt  conditions  (1.1)  in  a  simplified  setting  (to  avoid  nota- 
tlonal  difficulties).  Section  2  deals  with  preliminaries  and 
notation.  In  Section  3  we  formulate  an  abstract  approach  and 
based  on  it  prove  that  the  suggested  finite  element  formulation  of 
the  constraint  boundary  condition  leads  to  the  optimal  rate  of 
convergence  of  the  p-version.  Section  4  addresses  some  imple¬ 
mentations!  aspects. 
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2.  The  basic  notation  and  preliminaries. 


2.1.  The  Sobolev  spaces 
2 

Let  R  be  the  two  dimensional  Euclidean  space,  x  *  (x^.Xg) 

2  2 

e  R  .  Let  Q  c  R  be  a  bounded  Llpschltzlan  domain  with  the  boun¬ 
dary  r  «  dQ .  He  will  assume  that  r  Is  a  Jordan  curve,  r  * 
m 

U  f,  where  r.  are  smooth  open  arcs  with  parametric  description 
1-1  1  1 

Fi  *  Ux1*x2,|xl  “  xi(i<*)»  x2  "  xi , 2  ^  ^ '  l?l  <  1) ' 

1  —  1 , . . . ,m. 

Denoting  I  -  (-1,1),  r  ^  is  obviously  the  image  of  I  by  the  map¬ 
ping  Fj  -  (xA  1»*1  2>,  -  F1(I).  If  u ( s )  is  defined 

on  rA  then  by  D(?)  -  u(FA(?))  we  denote  its  transform  on  I. 

The  ends  of  1*^  .  will  be  called  vertices  and  denoted  by  - 

( xi , j ( — 1 ) ,  xlf2(“1,)'  Bi  *  ^Xi,l<1) ,Xi,2(1) ) *  Me  Wl11  further 

assume  that  B j  «  A1+1 ,  •  *  *  l,...,m.  By  this,  the  orien¬ 

tation  of  r i  Is  established.  In  general  we  will  denote  the  ver¬ 
tices  by  A^  (-  Bi_1),  1  -  l,...,m.  The  scheme  of  the  domain  and 
the  pertinent  notation  Is  shown  In  Figure  2.1. 


Figure  2.1.  Scheme  of  the  domain  and  notation. 
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Remark  2 . 1 .  We  assumed  that  the  domain  Q  is  simply  connected. 
This  assumption  has  been  made  only  for  notational  simplicity. 

Remark  2.2.  We  assumed  that  the  domain  is  Lipschitzian.  Once 
more,  our  results  are  valid  (with  proper  modification)  in  the  case 
when,  for  example,  some  arcs  coincide  (as  in  the  case  of  the  slit 
domain) . 

Remark  2.3.  We  have  assumed  that  the  arcs  are  sufficiently 

smooth.  For  the  sake  of  simplicity  we  assume  that  they  are  C°° 

arcs  (i.e.,  the  functions  x.  .,  i  =  l,...,m,  j  =  1,2  are 

1  f  J 

functions) . 

Jr 

By  H  (Q),  k  i  0  integer  we  denote  the  usual  Sobolev  space 

of  functions  with  square  integrable  derivatives  on  Q.  The  norm 

will  be  denoted  by  !j  •  j!  .  .  If  l  <  q  <  €+1,  l  >  0  integer, 

Hk(Q) 

then  we  define  Hq(Q  )  »  (H*(0  )  ,H<+1  (Q  )  )Q  ,  9  »  q~e  where  by 

( •  ,  •  )g  we  denote  the  usual  Interpolated  space  using  the  K-method 

(see  [5]).  The  scalar  product  (*,•)  and  the  norm  !  • :| 

Hq  Hq(Q ) 

are  defined  accordingly, 
k  — 

By  C  (Q ) ,  k  2  0  Integer,  we  denote  the  space  of  all  func¬ 
tions  with  k  continuous  derivatives  on  Q.  It  is  possible  to 
show  that  Hk(Q)  « — ►  C°(Q)  for  k  >  1,  where  by  < — ►  we  denote 
continuous  imbedding.  On  the  other  hand,  H^fQ)  <r  C°(0). 

b 

For  I«(-l,l),H(I),k^O  is  defined  analogously  as 
before.  If  k  >  1/2  then  Hk(I)  C°(i)  but  Hk(I)  <z  C°(i) 

for  k  s  1/2. 

V 

So  far  we  have  defined  H  (I),  k  i  0.  We  will  also  be  inter¬ 
ested  in  Hk(I),  k  <  0.  We  define  for  k  i  0 


[  uvdx 
J  -1 


H~k(I) 


veHk(I) 


V 

H*  ( I ) 


(Let  us  remark  that  sometimes  (see  e.g.  (5])  our  space  H  (I)  is 
denoted  by  (Hk(I))'  whereas  H~k(I)  is  used  to  denote  the  dual 
space  of  Hk  ( I ) ) . 

If  u  is  defined  on  then  we  define 

Hk(ri)  -  (u|u(Pi(J; ) )  *  U(?)  e  Hk  ( I )  > 

!!  ul!  „  -  !!U|i  v 

Hq(rA)  hk(D 

So  far  we  have  considered  only  scalar  functions  on  Q  and  I. 

The  spaces  of  vector  functions  are  defined  by  Cartesian  products, 

j  If  0  lr  2 

3H *(Q)  -  (HK(Q))  . 

Let  now 

Q  -  { ( Xj , Xj ) ( | Xj |  <  1,  |x2l  <  1) 

rl  *  <  < , ac2 )  |  | Xi |  <  1,  x2  -  -1} 

Q  will  be  called  the  standard  square  and  1  *  1,2, 3, 4  its 

sides  ( r i  ■  2,3,4  are  defined  analogously  to  r ^  in  an 
obvious  way) .  Let 

T  »  ((Xj.XjHIXjI  <  1,  0  <  x2  <  (1+Xj)  /5  for  Xj  <  0, 

0  <  x2  <  (l-Xj)/?  for  xx  >  0) 

T 

Y j  »  f(x  ,x  )!!x  !  <  1,  x2  *  o>. 

T  will  be  called  the  standard  triangle  and  y*,  i  ■  1,2,3  its 

sides . 


Let  us  remark  that  the  sides  of  T  and  Q  are  each  of  length 


2.  Later  we  Mill  often  not  distinguish  between  and  I. 


Figure  2.2.  The  scheme  of  the  standard  square 
and  standard  triangle. 

We  now  define 
2 

Pp(Q)  *  {u|u  is  a  polynomial  of  degree  s  p 

in  each  variable  and  x2  over  Q}. 

fp(T)  ■  {u|u  is  a  polynomial  of  (total)  degree  s  p  on  T) . 

Fp(I)  *  (u|u  is  a  polynomial  of  degree  s  p  on  I). 

We  have  then 

Lemma  2.1.  Let  v  e  n  (respectively  v  e  H*(?  *)  n 

flrh)  such  that 
P  1 


respectively 


i!v!!  _  s  p“(1_t)A,  t  =  0,1. 

IT(rJ) 

2  1 

Then  there  exists  u  e  P  (Q)  (respectively  9  (T))  such  that 

r  Sr 

u|  n  »  v  (respectively  u|  T  ■  v> ,  u|  n  =  0  (respectively 

r\  dQ-r J 

u|  T  «  0)  and 
dT 


!f u!:  1  S  Cp“1/2A 
H1  (Q) 

respectively 

1! ui;  1  s  Cp  1/2A. 
H1  ( T) 


For  the  proof  see  [2]  or  [3]. 


2.2.  The  model  problem 
Let 

2hJ(Q)  c  X(Q)  c  2H1(Q) 

where  tf(Q)  is  closed  in  aH1(Q).  JK(Q)  will  be  called  the  con¬ 
straint  space.  Assume  that  there  is  given  a  continuous  bilinear 

form  B(u, v)  on  2H1(Q)  *  2H1(Q),  u  *  (ui'u2^'  v  =  (vi'v2^  such 
that 


(2.1)  B(u,u)  *  ylluli2  .  ,  r  >  0  for  any  u  e  2H1(Q). 

H^Q) 

Then  obviously  for  any  Gj  e  (X(Q))',  there  is  a  unique  uQ  s  ^(Q) 
such  that 


B(uQ , v)  =  G^v) 

holds  for  any  v  e  )t(Q).  We  also  have 

'!  u  !!  s  C;!G  |! 

0  2H1(Q)  1  (2H1(Q)  )' 
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Denote  X  (Q )  *  (u  e  2 H1 (Q ) ,  u-p  e  X(Q)).  H  (Q )  will  be 

M  r 

called  the  p-hyperplane .  Then  our  model  problem  is  given  by: 

Find  u.  e  X  (Q)  such  that 

Op  . 

(2.2)  B(uQ,v)  =  G^v),  V  v  e  X  (Q  )  . 

We  have  then 


(2.3) 


II  u  'i 

°'2H1(Q) 


<  C(  !!p!! 


2H1(Q  } 


•IG  11 

YV 


(Q) )' 


If  p  =  0  then  we  will  speak  about  a  homogeneous  constraint 
problem  while  for  p  *  0  we  will  speak  about  a  nonhomogeneous 
constraint  problem.  We  call  these  constraint  problems  because 
X (Q)  *  2H1(Q) . 

There  are  many  constraint  problems  in  applications.  We  will 
consider  the  one  when 


(Q)  *  <(u1#u2)  e  2H1(Q)|  ujr  *  0,  k  =  1,2,  j  =  1 . m) 

J 


-e=i 


3 


(3) 


(ak  are  matr*ces  smooth  functions  on  (say 

( 3 ) 


where  a 

^(fj)).  Additional  assumptions  on  (a^J^)  will  be  imposed 
later. 

Obviously  when  k  =  1 '  ak  €  “  0  for  ^  *  *  we  9et  Diri_ 

chlet  boundary  conditions  (in  general  we  get  Dirichlet  conditions 

3 


have  rank  2  for  all  x  e  f  ) 


when  a 

is  no  constraint  and  we  have  the  Neumann  problem. 


If  (a^]}  =  0  there 


,3 


If  aJ  has  rank  1  then  we  can  write  the  constraint  on 


as 


cY>  u  +  a(^  u 

1,1  1  1,2  2 


0 


which  will  be  written  in  the  form 


(2.4a) 


a  '  “"u1  +  ;?  '  J  'U2  =  0  . 


Obviously  if  P  =  iP1,P2)  then  the  nonhomogeneous  constraint 
problem  is  characterized  on  f j  by 


(2.4b) 


a^u,  +  *<J>u2  =  a(i>Pl  +  3'»p2. 


Problems  of  this  type  are  common,  for  example,  in  the  theory 
of  elasticity.  For  simplicity  of  the  exposition  and  notation  we 
will  restrict  ourselves  to  the  model  problem  where 


B(u,v)  = 


(2.5) 


dv1  dn1  dv1  <3u2  dv2  dn2  dv2 
flx'j  c>x  +  dx2  dx2  +  dx1  3x^  +  dx"  0x2 


+  ulvi  +  U2V2  dxidX2 


G^v)  =  ( f i  +  f2v2)dx1dx2,  f  =  ( f a , f 2 )  c  H  (0). 


Although  we  restrict  ourselves  to  this  special  case,  our  results 
hold  in  general,  e.g.,  for  elasticity  problems,  etc. 

(  1  )  X  — ■ 

We  will  assume  that  ct^  ^  ‘2  C  (Fj).  In  practice  we  have  the 
nonhomogenous  constraint  problem  defined  so  that 

i)  If  )  has  rank  2  on  then  the  constraint  is 


Zc,k;?J  U<  =  9k’‘  k  * 


where  (g|^ ,g2^ )  are  defined  on  f  ^ .  Hence  obviously  we  can 
transform  the  above  constraint  equation  to 


r*»: 

|*i 

l** 

i*i 

t*> 

[r 


1 


Because  are  assumed  to  be  smooth,  we  see  that  (g^^),  i  = 

1,2  have  the  same  smoothness  as  (gj^). 

ii)  If  (et^  j  }  has  rank  1  then  the  constraint  equation  is 

=  g(^. 

We  add  the  condition 

a ( j )  +  j ( J )  0  o  on  f  . 

This  enables  us  to  transform  the  constraint  equation  to 

(2.7)  a ^ ^ u.  +  3 ( 3) Ug  =  g ^  ^  with  u  ^  ^  +  j ( J )  -  i # 

As  formulated  above,  g ^  ^  are  defined  separately  on  each 
fj.  We  will  assume  that  g^^  satisfy  consistency  conditions, 
namely  that  there  exists  p  =  (P1,P2)  s  2H1(Q)  such  that 

(2.8a)  ^  =  p  ja-  ,  =  p  |p  , 

l  1  1  j  2  2  j 

respectively 

(2.8b)  Ci(j)  p  |r  +  ?(J5p  |p  =g{j). 

j  j 

These  conditions  have  to  be  imposed  especially  at  the  vertices  of 
0  . 

The  sides  where  the  constraint  (2.6)  is  imposed  will  be 

called  total  constraint  sides,  while  where  the  constraint 

(2.7)  is  imposed  will  be  called  partial  constraint  sides:  We  will 

enumerate  the  total  constraint  sides  as  f  .  ,  j  =  l,...,m  and 

j  1 

the  partial  constraint  sides  as  f  ,  j  =  m  +l,...,m. 

i  1 

2.3.  The  p-verslon  of  the  finite  element  method 

Assume  that  the  domain  Q  has  been  partitioned  into  a  finite 
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n 

number  of  subdomains  0.,  i.e.,  Q  -  U  0  ,  We  shall  assume  that 

1  1=1  1 

Qj  is  the  curvilinear  quadrilateral 


or  curvilinear  triangle 


Q 


7i  (Q) 


Q 


' ri(T ) 


where  Q  and  T  are  the  standard  square  and  triangle,  respective¬ 
ly.  The  domains  Q ^  will  be  called  elements.  We  will  assume  that 
f  ^  is  a  smooth  one  to  one  mapping  of  0^  onto  Q,  respectively 
T.  It  is  obvious  what  the  vertices  and  sides  of  Q ,  correspond 
to.  If  /  is  a  side  of  then  induces  mapping  of  I 

onto  ;  (realizing  that  all  the  sides  of  the  standard  square  and 
triangle  have  the  same  length  as  I). 

We  shall  assume  the  following  about  the  partition  and  tfle 
mappings  : 


.  a.)  If 

°i  °J 

ii 

M- 

3'  0 

vertex  or  a 

side  of  both 

0  ^  and 

3 )  If 

Ri.j  *  > 

i<3 

then  we 

of  I  onto 

/  ^  .  induced  by  the 

tical.  We  denote  F 

by 

Fi.r 

A,B  the  vertices-  of 

Qi 

and  0  j 

Assume  that 

(a1.b1) 

and 

( a2 ' b2 

Q  T 

•vk  or 

such  that 

Then  if  C  e 

r  ^  ,  and  C 

-  »1(c1 

= 

is  either  a  common 


Since  we  assumed  that  are  smooth  mappings,  the  vertices 

of  Q  necessarily  have  to  coincide  with  some  of  the  vertices  of 
Q^.  We  will  further  assume  that  for  any  Fj  there  is  an  element 
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Q1  such  that  one  of  its  sides  coincides  with  r  ^ .  This  assump¬ 
tion  is  made  without  any  loss  of  generality. 

Denote  now 

i'PlQ)  •  {u  -  H 1  ( Q )  u|  (O)  «  ^(Q)  if  is  a  q^adri- 

P  u  i  1  p  1 

lateral  and  u|Q  (T  i(?)  )  --  f*(T)  if  Q  i  is  a  triangle}, 
**  (0)  *  { 1  tPp  (Q  )  )  2  . 

2P  (/)  -  {u  u.|,(F.(t))  €  f>  (I),  j  -  1,2}. 

P  J  /  1  P 


Here  }  =  T  ^  or  any  side  of  an  element.  Let  us  define  the  con¬ 
straint  space  Vp(0 )  c  2,Pp( 0)  as  follows. 


i) 

If 

is  a  total  constraint  side  with  end 

points 

Aj'Aj+l 

and 

u  s  X9^(Q),  then  *  ui(Aj  +  i)  * 

o 

>-*• 

H 

*-» 

and 

u.c.ds  ■  0,  for  all  n.  €  lf>  _ ( r  ) ,  i  *  1,2. 
i  i  i  p-2  J 


ii)  If  Tj  is  a  partial  constraint  side  then 
(o(J)u1  +  l(3)u2)(Ak)  -  0,  k  -  j.J+1 

and 

(a^^u1  +  ?^^u2)tds  ■  0,  for  all  «r  s  2Pp_2(r  j)  • 

K 

The  p  hyperplane  *  (Q)  is  defined  analogously.  Let  g 

r  P 

be  defined  in  terms  of  p  by  (2.8).  Then  on  we  impose 

ui(Ak)  *  gfJ)(Ak)'  1  *  x’2*  k  “  ;)'j  +  1 

respectively 


and 


(o<j)u1  ♦  '?(J,u2)(Ak)  -  g,J)|Ak),  k  -  j.j+1 


(2.9a) 

raapactivaly 

(2.9b) 


ui*ids 

Jr , 


gJJUida,  i  ■  1,2 


J 


(a  *  ^Uj  +  H  *  ^u2  )?ds 


g(*W 


J 


J 


The  p-varsion  la  than  daflnad  analogoualy  aa  bafore. 
(Q)  auch  that 

P  p 


Find  u. 


(2.10) 


B(u  ,v)  -  G.(v)  V  v  e  >  (Q) . 
pi  p 


Raaark .  Constraint*  of  tha  typa  consldarad  ara  typical  In  elasti¬ 
city  thaory.  Hara  and  u2  ara  tha  dlsplacaaants  In  the 

directions  Xj  and  x2  raapactivaly .  As suae  now  that  the  dis- 
placeaent  Is  constrained  In  tha  noraal  direction  only  (and  Is 
friction  free  In  tha  tangential  direction).  Than  on  the  boundary, 
we  obtain  tha  partial  constraint  Uj  cos  a  ♦  u2  sin  a  »  0  where 


a  la  tha  angle  of  tha  outer  noraal  with  tha  axis  x^ . 


3 .  The  convergence  of  the  p-vtr«lon  of  the  finite  element  method 


3.1.  An  result 

Mt  will  first  dsscrlbs  an  abstract  framework  which  will  be 

the  basis  for  the  forthcoming  analysis. 

Let  X  and  H  be  Hilbert  spaces  and  X  c  x,  w  c  w,  p  • 

P  P 

1.2....  be  one  parameter  families  of  finite-dimensional  subspaces 
Xp  c  x,  Wp  c  w  will  denote  corresponding  families  of  hyperplanes 
such  that  (u-v)  «  X_  whenever  u.v  e  X„  and  («-*>)  e  w  when- 


Let  a(u,v).  u.v  c  x  be  a  continuous  bilinear  fora  on  XX 
and  b(v.e)  be  a  continuous  bilinear  form  on  X*W  such  that 
b(v,e)  S  Cl!vi>x!‘*:!w. 

Let  u_  «  X,  c  W  and  u_  *  X  .  •  «  W_  be  such  that 

0  0  P  P  P  P 

(3.1a)  a ( uQ , v )  ♦  b(v,*Q)  -  F^v)  V  v  e  Xp 

(3.ib)  b(u0,e)  ■  f2(r)  v  r  «  wp 


and 

(3.2a) 

(3.2b) 

Define 

have 


a ( u  . v )  ♦  b(v,e  )  -  F  (v)  V  v  c  X 
P  Pi  P 

b ( u  , * )  -  F_ { ? )  V  m 

p  2  p 

Z  -  {vc  X,  b(v,a)  ■  0  V  *  e  W  )  c  x  Then  we 
P  P  P  P 


Theorem  3.1.  Let  a(u,u)  «  filu11*  ,  r  >  0,  for  any.  u  e  Z  .  Then 

x  p 


(3.3)  !i UQ~up  x  s 


inf 

b( un-w  ,  * )  *0  Ve*S 

op  P 

u  -w  <X 
P  P  P 

(i.e.,  u  -w  «Z  ) 

P  P  P 


u_-w  „  ♦ 
0  p  X 


inf 

t  ew 

P  P 


*0_lp  W1 


15 


Proof,  for  arbitrary  *»  «  X  and  r  €  W  wa  have 

P  P  P  P 

(3.4a)  a(u  -w  .v)  ♦  b(v,e  -t  )  ■  a(urt-w  ,v)  ♦  b(v,*.-»  )  V  v  e  x 
PP  PP  °P  °P  P 


(3.4b) 


b(u  -w  ,e)  »  b(un-w  .*)  V  r  e  W 
P  P  Op  P 


For  v  •  Zp  (3.4a)  ylalds 

(3.3)  a(u  -w  , v)  -  a(u  _-w  ,v)  ♦  b(v,* n-r  ) 

P  P  Op  Op 


Suppose  now  that  w  Is  such  that 

P 


Then  by  (3.4b) 


and  hence 


b(Wp.r)  -  b(uQ,r)  v  *  e  wp. 


b(u  -w  ,*)  -  o  v  *>  c  w 
p  P  p 


u  -  w  e  Z. 
P  P  P 


Mow  using  v  ■  up  -  wp  In  (3.3)  we  get 

•<VVW  !  Ct!oo-Vx  'W'x  *  »o-*p'h  V"p  X1 

and  hence  by  coerclvlty  of  a(*,*)  on  Zp, 


it u  I*  $  Cf '  u  -m  u  ♦  ii 0  —  r  1 
P  p‘ X  1  0  p’ X  (p  NJ 


and  hence  also 


::Ua-w.L  s  C[  urt-w  •'  ♦  '!S«-T  I-.] 

U  P  A  v  p  A  U  p  “ 

fros  which  (3.3)  follows. 

3.2.  The  convergence  of_the  p- version 

Lat  uQ  »  (uQ  j,uQ  2)  *  aHl(Q)  be  the  solution  of  our  con¬ 
strained  problea  (2.2),  (2.3)  and  u  »  (u  ,,u  .)  e  2?  <Q)  be 

P  P« 1  P> 2  p 

the  approxiaatlon  given  by  (2.10). 
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flu. 


tin  will  mium  that  uQ  c  aH*(Q ) ,  k  >  3/2.  Hanca  y--  « 
aM°(r1) .  i  ■  1,2,..., a.  Lat  tha  constraint*  on  Tj  ba  as  in  sac- 
tlon  2.2,  with  r.  ,  J  -  1 . a.  balng  tha  total  constraint  sidas 

J 

and  r.  ,  J  *  a1’tl....,a  tha  partial  constraint  sidas.  Than  It 
can  ba  varlfiad  that  for  any  v  <  aH1(0>, 

"‘l 


®(u, 


SViv  .  !Vlv' 
an  V1  an  v2 


J-1  'r 


ds 


(3.«) 


[  (a  ♦  *  Jv2) 

rduo.i 

9H 

an 

an 

L  J 

ds  •  G j ( v ) 


J 


whara  wa  hava  asauaad  a  J  «■  n  J  -  i.  Moraovar,  for  s.Sj.Sj  and 
p  aa  In  (2.6)  -  (2.9), 


Z  [  ‘vi'i  *  up.j’2,d*  *  Z  I  (al\.i  *  ,ljup. 

3-i  Jr  J-,*‘  Jr. 


)  eds 


(3.7)  -  £  ,u0.l'l,'U0,3,3|d*  *  Z 


(o  Juo.l*i  u0,2l*<i 


j-i  r 


J-«,*i  Jr 


*1  ,  ■  r  1  1 

y  (plpl','p2f2)d#  *  X  (a  ipl^  *f>2)9du 

J-1  Jr  j-^1  Jr 


and 


(3.6a)  up  k(Aj)  ■  uQ  k(Aj)  •  Pk(Aj),  J  -  «.*♦!.  k 


-  1.2, 


(3.6b)  *  (a'pjKA.)  *  (i/pj)(A,) 

)  -  /./♦!.  <  -  i.i+1 . 1.. 

( Wa  rturk  that  uQ  ^(Aj),  i  ■  1,2  has  meaning  bacausa  wa  assumed 
that  uQ  «  3 1^(0) ,  k  >  3/2)  . 

Ns  now  dsflna 

X  -  aH1(0)  ,  «•«  -  !!•!!_  . 

A  aH  (0) 

and  for  any  6  ■  (Sj.Sj)  €  X, 

XP,6  *  <u  “  (VU2}  c  2%(Q)*  W  “  6k(V'  k  "  l'2 

J  -  i  -  lj . i  .  (a^UjMAj)  ♦  (ft*u2)(Aj) 

*  J  -  /  »  i  +1 . i,,). 

Wa  than  taka  in  our  abstract  framawork 

X  ■  X  ,  X  -  X  _ 

P  p.p  p  P.0 

whara  p  satlsfias  (2.6)  -  (2.8).  Moraovar,  1st 

W  -  n  aH_1/2(r  )  .  n  lH'1/2(r  ) 

j-1  3-a^i 

with  tha  norm 


Wa  saa  than  that  aQ  *  W  whara  (saa  (3.6)) 
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.  (H.n  aV 
>,3  9n  L  '  5n 

1  *j 

.  faVi  .  au0.2]| 

~Sn~  *  -3n~  I 

to  * 


i  3  *  1 1  • .  •  i®< 


Define  W  -  W  «  n  c  w 

P  P  jUiJ  P 


,  3  -  *.+1 . m. 

Js 


where 


jWp  -  aH-1/2(r1^)  O  2*p_2(r3)*  3  -  1 . 


-  1H‘1/2(T1  )  «  J  “  ®i  +  1 . . 

3 

Let  •(•,*)  end  b(*,*)  be  bilinear  fores  defined  respec¬ 


tively  on  X*X  and  X*W  by 


a(u, v)  -  B(u, v) 


"I-  ■  r  i  i 

i.P)  ■  £  ♦  u2a2)do  ♦  (a  +  9  ^u2 )  e 

3-i  Jr  3-^1  r. 


It  may  be  seen  that  the  right  hand  side  of  (3.7)  defines  a  linear 
functional  G.  on  W  .  Then  (3.6)  -  (3.8)  show  that  (u_,e-) 

2  p  0  0 

satisfy  (3.1)  with  ■  G^,  k  ■  1,2.  Moreover,  if  we  can  find  a 
unique  pair  (up,#p)  satisfying  (3.2),  then  u^  will  be  precise¬ 
ly  our  finite  element  solution  satisfying  (2.10).  He  will  now 
verify  that  the  mixed  method  defined  above  satisfies  the  assump¬ 
tions  of  Theorem  3.1.  This  in  turn  will  lead  to  the  existence  and 

uniqueness  of  the  solution  (u  ,*  )  of  (3.2)  and  an  estimate  of 

P  P 

the  rate  of  convergence  of  up  to  u. 

Obviously,  a(u,v)  satisfies  the  desired  continuity  and 


R: 


•t 

Ev 


•?T 
t  *• 


Kf 

$ 

S 


:wn 


coerclvity  conditions.  For  J  ■  1,2 


■  ws  have 


u r  ds 


'j 


^  u  1/2  -1/2 

HA/,i(r  )  H  *  (V  ,  ) 

*1  i 


from  which  the  continuity  of  b(*,*)  may  be  deduced.  Hence 
Theorem  3.1  is  applicable.  Let  us  now  estimate  inf  <p  -  t 


rsW 


W 


P 

.k-3  /  2 


First,  let  m.,  +  1  s  j  i  m.  We  assumed  that  o  .  :  H~  '"“(T.  ),  k 

l  J  1 j 

lc-3  /  2  i  (  T  ' 

3/2.  Hence  oQ  (FA  (<  )  )  *  *(')  «  HK  J/<£(I).  Let  a  -z  'fp_2 

J 

be  such  that 


(3.9) 


op  d* 


*p  d /  ,  V  p  e  YjI!)  . 


Then,  with  q  »  *  -  a,  we  have 

- (k-3/2 ) 


(3.10) 


HW(I) 


s  Cp 


Hk-3/2(i) 


Now,  for  arbitrary  v  €:  H  (I),  we  have  by  (3.9) 


}  .qv  d?  j  qfv-a^df  q 


_  v-a ,  _ 

I  *  H°  (  I )  1  ”° 


H^I) 


H1  ( I ) 


v 


Oil  £  Cp'1  q 


H1  ( I ) 


H°  ( I ) 


where  <r  is  a  polynomial  of  degree  p  -  2  satisfying 

v  -  ai  o  -  Cp_1  v  1 

A  H  (I)  H  (I) 


This  yields 
(3.11) 


-i 

H  1  ( I ) 


cp-(k'1/2)* 


Hk  3/2 (I) 


Interpolating  (3.10),  (3.11)  and  using  the  fact  that  F  is  a 

J 

smooth  mapping,  we  obtain 


inf  U, 


*€  W 

P  J  P 


V  -1/2  *  Cp 

p  h  '  (r .  ) 


-(k-i) 


°,j'Hk-3/2(r  ) 

J 


We  get  similar  estimates  for  r.  ,  j  =  l,...,m1>  so  that 

j 


(3.12) 


inf  ||0  -  .t  !!  s  Cp“(k“1)l!u  I! 

X pe Wp  0  P  w  0  2  n  (Q ) 


We  now  estimate  inf  !u  -  w  ||v.  Using  the  results  from  [2],  there 

O  P  A 


exist  z,  e  9  (Q )  ,  i  ■  1,2  such  that 
i  P 


(3.13a)  !!  u_  .-zJ|  s  Cp" ( k_t } !!  u.  >!  ,  i  =  1,2,  t  =  0,1 

0,1  1  HT(Q)  0,1  HK(Q) 


(3.13b)  uQ  a(N)  =  zi(N)  for  each  node  N  of  the  mesh 


(3.13c)  liu.  . -z . !!  s  Cp“(k"1/2-t)!!u_  J!  .  ,  t  =  0,1,  i  =  1,2 

1  Ht(r^)  0,1  HK(Q) 


j  *  1 , 2 , . . . ,m 


j  2 1 

Let  n^+1  s  j  s  m.  Let  us  denote  *  *  a  (uQ  j-Zj)  +  !i  (uQ  2~z2^  ’ 


Then  we  have  *(A.  )  *  *( A.  )  =0.  Let  x(Z)  *  *{F(t;))  and  let 

J  j 


t(F(S))  -  T(JJ)  e  **^(1)  satisfy  t(£1)  *  0  and 


*  * 

TV  cK  «  ?V  d(  ,  V  V  6  lF  2(I) 

"  T  T 


Because  of  (3.13b)  we  can  write 


r\y  d*  ■  x'u)' 

I  I 


di;  ,  V  m  e  lfp(I)  ,  w(±  1)  *  0 


and  hence  by  Lemma  3.2  of  [2] 


Hsr  -  T!i  s  Cp"(k"1/2-t,i!u  I!  .  ,  t  =  0,1 

IT(I)  U  2HK(Q) 


Using  (3.13c),  this  gives 


F 

k 


& 

B 


I-# 

8 

: 

i 


I 


i** 

[»*, 

[»* 

ft* 


is' 


.  *  *  »  i.i  r.  * 


- T •'!  t  s  Cp"(k~1/2"t,!iU  ji  .  t  *  0,1. 

H  (I)  H  (Q) 

Now  using  Lemma  2.1  it  follows  that  there  is  a  w  e  1f>  (Q)  such 

P 

that  w  *  0  on  Q  -  Q  where  Q  is  the  element  with  the  side 

j 


r .  ,  w  *  r  on  r.  ,  w  *  0  on  dQ-r.  and 


(3.14) 


!lw!!  *  Cp“(k"1>|!u  *! 

Ha(Q)  u  2HK(Q) 


Letting  w.  =  (w,w)  e  2H1(Q),  we  see  that  w.  will  satis- 

i 

fy  (3.14)  with  H1(Q)  replaced  by  2H1(Q). 

Let  now  w  =  z  +  wA  =  +  (w,w).  Using  (3.13b)  and 


J 


the  fact  that  w(A^)  =0,  we  obtain 
*  j  *  h 

(c  J(u0|1-wp  l)  +  fl  J(u0,2-Mp,2))(A<>  -  °-  e  *  VV1 


Moreover,  for  v  e  V  (T  .  ) 

p—2  1  j 


’3 


*  j 

(a  J(un  ,-w^  ,)  +  J(un  0-w_  „))?  ds 


0,1  p,l 


0,2  p , 2 ‘ 


J 


*  j  *  j 

(#  -  (a  J  +  (i  J)w)p  ds  « 


(#-w)f  ds  *  0 


j 


j 


*  j  *  j 

where  we  have  used  aJ+/3J  =  l.  We  may  construct  w,  as 

above  for  all  partial  constraint  sides.  An  analogous  construction 
can  be  carried  out  for  total  constraint  sides  as  well.  Then  if 

m 

i 

w^  ■  z  + 

P 

we  see  that 
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and 


u  -  w  e  X 
P  P  P 


b{u0-Wp,f/>)  =  0,  V  ip  e  Wp 


m 


|U0‘Wp:!X  '  liuo"z!lx 


+  Z!lwi»jt  s 

J- 1  3 


„  - ( k-1 ) 

<  Cp  !!U 


0  2Hk(Q) 


This  provides  a  bound  for  the  first  term  in  the  right  hand  side  of 
(3.3) .  Hence  we  have  proven 

Theorem  3.2.  Let  uQ  e  2Hk(Q),  k  >  3/2.  Then 


0  P,,aHi(Q) 


_  -(k-1),,  „ 

5  °P  ilUo!i2  k 

H  (Q) 


where  is  the  exact  solution  and  u  is  the  finite  element 

0  p 

solution  of  the  constrained  problem,  provided  that  uQ  and  up 
exist  .o 


The  next  theorem  deals  with  the  question  of  existence  and 
uniqueness  of  (uo,0o>  and  (up,0p) . 

Theorem  3.3.  The  (exact)  solution  (uo,<po^  of  the  constrained 

problem  exists.  The  finite  element  solution  (u  ,<p  )  exists  and 

P  P 

is  unique. 

Proof .  In  Section  2.2  we  have  shown  that  uQ  exists  and  hence 

(uQ,0o)  exists,  too.  The  finite  element  solution  (up'^p^  is 

determined  by  the  solution  of  a  linear  system  of  equations  with 

square  matrix.  Hence  the  existence  follows  from  the  uniqueness. 

Assume  therefore  that  there  is  a  solution  (u  ,v>  )  of  the  trivial 

P  P 

problem.  Obviously  u  **  0  *  0  is  also  a  solution  of  this  problem. 

Hence  u  *  0  because  of  Theorem  3.1.  We  have  to  show  therefore 
P 

that 
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(av:  +  d?  =  0 


implies  v>p  *  0.  Because  o  +  ,3  -  1  we  also  have  J  v</> ^  d?  =  0 

for  all  v  e  lf  (I)  r>  H*(I)  while  <p  e  x9  _(I).  This  leads  to 

p  0  p  p-2 

<*>p  =  0  which  leads  to  the  desired  result.  : 

Remark .  We  have  dealt  only  with  a  model  problem.  It  is  obvious 

that  the  theorem  holds  in  general,  as  for  example,  for  the  theory 

of  elasticity. 


Some  aspects  of  implementation. 


Here  we  will  make  some  comments  about  the  implementation  in 

* 

the  framework  of  the  code  PROBE  (see  [6]).  The  shape  functions 
are  defined  as  usual  on  the  standard  square  or  triangle.  There  are 
three  types: 

a)  the  model  shape  functions  which  are  linear  on  every  side 
of  Q,  respectively  T; 

b)  the  side  functions  which  are  zero  at  the  vertices  of  Q, 
respectively  T  and  on  y  are  of  the  form 


* 


j 


X 

)dx,  j  =  1,2, .  . . 

-1 


where  is  the  Legendre  polymonial  of  degree  j. 

is  then  a  polynomial  of  degree  j+1 ; 

c)  The  internal  shape  functions  which  are  zero  on  dQ 
(respectively  dT) . 

The  stiffness  matrices  are  first  computed  in  the  standard  way 
without  constraints.  Then  the  constraints -are  imposed  at  the  ver¬ 
tices  A..  This  only  involves  the  amplitudes  for  the  nodal  shape 
J 

functions.  Then  the  conditions  (2.9  a,b)  only  involve  amplitudes 
for  the  side  shape  functions.  The  functions  tp  in  (2.9  a,b)  are 
computed  as  derivatives  of  the  Legendre  polynomials  from  the  usual 
recurrence  formula  and  the  integration  is  made  using  numerical 
quadrature . 

The  condition  (2.9a)  is  especially  simple  because  u^,  = 
£c^j.  Integrating  by  parts  and  exploiting  orthogonality  of  the 


* 

The  code  PROBE  is  the  code  of  Noetic  Tech.,  St.  Louis. 
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rn»  MHCttn  HC.  JMeiSi  JOUXHiL  *•  m  t»Ur»l  part  of  tba 
XnllMi1  for  n^ttiuloliMi  aad  Tinhnalogy  of  febo  University  of  Maryland, 
under  the  paaeral  adalal  strati  on  of  Um  Diroetor ,  Institute  for  Physlaal 
Solanoa  and  Taohaology.  It  boo  tbo  following  900101 

o  To  conduct  manor  sh  la  tbo  aatbaaatloal  theory  and  oaaputatlonal 

lapl adaptation  of  auMrloal  analyala  and  related  top  lea.  with  eophaels 
oa  tba  auaorloal  treataoat  of  llnaar  and  nonllnoar  dlffarontlal  equa¬ 
tions  and  problana  la  llnaar  and  nonllnoar  algebra. 

o  To  balp  bridge  gaps  between  oaaputatlonal  dlrootlona  In  aaglnaorinc, 
physloa,  a  to.,  and  those  In  tba  aatbaaatloal  oo—lty. 

o  To  provide  a  Hal  tad  nonsuiting  aervloe  in  all  araaa  of  nuaarloal 
aatbaaatloo  to  tba  Uhl  varsity  an  a  whole,  and  alao  to  govornaant 
agoaolao  and  laduatrlaa  In  tba  Stata  of  Maryland  and  tba  Washington 
Matropolltan  araa. 

o  To  aaalat  with  tba  aduoatlon  of  nuaarloal  analyata,  aapaolally  at  tha 
poatdootoral  level,  In  ooajuaotloa  with  tba  Zntardlaclpllnary  Applied 
Math  Motion  Prograa  and  tba  progroaa  of  tba  Hathoaatloa  and  Coaputar 
Solanoa  Papartaanta.  this  lneludoa  notion  oollaboratlon  with  go vern- 
aaat  agoaolao  aueh  aa  tba  Motional  Burnau  of  Standarda. 

o  To  bo  aa  intaraatlonal  ooatar  of  atudy  and  monarch  for  foralgn 

atudanta  in  warloal  aathonatloa  da  ara  aupportad  by  foralgn  govern- 
aoata  or  exohaage  aganolao  (Pulbrl^t,  oto. ) 

furthar  lnforaalloft  any  ba  obtalnad  tram  Profoaaor  X.  Babulka,  Chalman, 
Laboratory  for  Nuaarloal  Analyala,  Institute  for  Physical  Solanoa  and 
Taohnology,  University  of  Maryland,  Collage  Park,  Maryland  207*2. 
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